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Abstract. - The methods of statistical mechanics are applied to two-dimensional foams under 
macroscopic agitation. A new variable - the total cell curvature - is introduced, which plays 
the role of energy in conventional statistical thermodynamics. The probability distribution of 
the number of sides for a cell of given area is derived. This expression allows to correlate the 
distribution of sides ("topological disorder") to the distribution of sizes ("geometrical disorder") 
in a foam. The model predictions agree well with available experimental data. 
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Introduction. — Foams and related physical systems 
(like emulsions, biological tissues, or polycrystals) are 
ubiquitous, and serve as a paradigm for a wide range 
of physical phenomena and mathematical problems [IH1]. 
One of them deals with the general topological and ge- 
ometrical properties of cellular materials [5-12 j. In this 
connection, Quilliet et al. [5,6 studied recently the topo- 
logical features of two-dimensional (2D) soap froths un- 
der slow oscillatory shear. Such macroscopic strain in- 
duces rearrangements within the foam, and the number of 
sides of every cell evolves in time through local topological 
changes (Tl events) [Tj [T3l [Tl]. Nevertheless, Quilliet et al. 
reported the existence of an equilibrium state after few 
cycles, characterized by a stationary probability distribu- 
tion of the number of sides per cell (topological disorder). 
They also showed that the width of this distribution of 
sides is strongly correlated to the distribution of bubble 
sizes (geometrical disorder) within the foam. These results 
suggest that the macroscopic state of a homogeneously 
sheared foam can be adequately described using the ideas 
and formalism of statistical thermodynamics. Indeed, the 
pioneering work of Edwards on granular matter [15] has 
shown how the powerful arsenal of statistical physics can 
be extended to athermal systems. This method has proven 
its applicability to other fields as well [IB] . Because there 
is no thermal averaging due to Brownian motion, this ap- 
proach requires the presence of a macroscopic agitation 
(analogue to an effective temperature) allowing the sys- 
tem to explore its entire phase space. 



Various attempts have been made in the past to de- 
scribe the geometrical and topological properties of 2D 
foams using the concepts of statistical thermodynamics 
[TTH22] . However, these former theoretical approaches rely 
on strong assumptions: either they use an ad- hoc inter- 
action potential between bubbles [UJ[T8], involve (rather 
than deduce) empirical laws correlating size and side dis- 
tributions [T9H2T] , or ignore some geometrical constraints 
(for instance, only the mean bubble area is specified, not 
the individual bubble areas [19,20,22 ). Some of these 
models are based on the maximum entropy (information 
theory) formalism [T9] |20], which has been subject to con- 
troversy [23]. Other models invoke minimisation of the 
energy [21], or a combination of both principles [22] to 
describe the state of a foam. However, it has been estab- 
lished 12.21 127] that different arrangements (topologies) 
of a large number of bubbles do not really affect the energy 
(see discussion below). Furthermore, none of these models 
can account for the correlations between topological and 
geometrical disorders reported by Quilliet et al. [5j[6]. 

In this letter we set up a framework for describing the 
equilibrium state of a two-dimensional foam, basing our 
development on analogies with conventional statistical me- 
chanics. As for other athermal systems [15J[T6], we show 
that the energy is not relevant to describe the macroscopic 
state of a foam. Instead, a more appropriate state vari- 
able is introduced: the total cell curvature. We establish 
the function of state which is minimized for a finite cluster 
of bubbles at equilibrium. This thermodynamic potential 
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function differs from entropy or energy used in previous 
theories. The formalism developed here allows to derive 
an analytical expression for the distribution of the num- 
ber of sides. We show that the semi-empirical expression 
conjectured by other authors [17, 18 , 28 , 29] is recovered, 
in a certain limit. Finally, the bubble size-topology corre- 
lations deduced from the present theory are investigated, 
and compared with experimental data. 

Physical, geometrical, and topological con- 
straints. — Consider a given set of Nb bubbles with pre- 
scribed areas {Ai} (we focus on time scales much shorter 
than those typical of bubble coarsening and coalescence, 
so the bubbles preserve their integrity and size 1 ). A 2D 
foam is a partition of the plane without gaps or overlaps, 
and its structure must obey certain constraints pQ (one 
considers the dry foam limit where liquid volume fraction 
is negligible). The physical constraints follow from the 
mechanical equilibrium throughout the system: first, the 
balance of film tensions at every vertex implies that the 
edges, or Plateau borders, are three-connected making an- 
gles of 120° with each other (Plateau's laws). Then, the 
balance of gas pressures in adjacent bubbles implies that 
every edge is an arc of circle, whose algebraic curvature 
Kij is proportional to the pressure difference between the 
two adjacent bubbles i and j (Laplace's law): 

Pj — P% 

Kij = — Kji = - , (1) 

where 7 is the film tension, and Pi and Pj are the pressures 
in bubble i and j, respectively (by convention, Kij > 
when the center of curvature is outside the cell z, i.e.: when 
Pj > Pi)- As a consequence, the algebraic curvatures of 
the three edges that meet at the same vertex must add to 
zero: 

= 0. (2) 

The foam must also satisfy topological and geometri- 
cal requirements: apart from the constraint of prescribed 
areas, its structure must obey Euler's rule, which relates 
the number of bubbles Nb, with those of Plateau borders 
7Vp5, and vertices N v : Nb + N v — Np^ = c, where c is 
the Euler-Poincare characteristic (c = for a torus, c = 2 
for an infinite plane). This rule, combined with Plateau's 
laws, immediately gives: 

N Pb = 3(N B - c). (3) 

Finally, the Gauss-Bonnet theorem applied to a n-sided 
cell yields Y^j=\ kj^ij — ^{n — 6)/3, where Uj denotes the 
length of the edge common to bubbles i and j. 

Microcanonical ensemble (iV#, Ktot, N 3 ). — In sta- 
tistical mechanics, the most fundamental entry is certainly 

1 Rigorously, the area and pressure of a cell vary from one config- 
uration to another; only the number of molecules of gas it contains 
remains unchanged. We consider that the area fluctuations are small 
so that each bubble can be identified by its area. 



via the microcanonical ensemble. Suppose that the foam 
is agitated slowly as compared to the characteristic relax- 
ation time after a Tl process [T|[T3|[T^] (one considers only 
perturbations that preserve the area of every bubble, and 
so the total foam area). Then, the deformation of the foam 
is quasistatic: its structure evolves through configurations 
which always satisfy the physical, geometrical and topo- 
logical constraints stated above. Such configurations will 
be referred to as accessible micro states. 

By analogy with the fundamental postulate of statisti- 
cal mechanics [30], we hypothesize that all the accessible 
microstates of a given set of bubbles filling the 2D space 
have equal probability, under a slow macroscopic agita- 
tion. Obviously, a periodic or infinite foam fills the 2D 
space. It must be pointed out that the postulate applies 
to a free cluster too, provided that the surrounding air is 
included as a supplementary bubble. Then, this "extra 
bubble" must also be taken into account in the counting 
of the number of accessible microstates. In the following, 
an infinite or periodic foam, or a free cluster plus the sur- 
rounding air shall be referred to as an unbounded foam. 
Conversely, a free finite cluster, a cluster within a larger 
foam, or a foam enclosed in a container shall be referred 
to as a bounded foam. Figure [1] summarizes the different 
kinds of boundaries that exist for a 2D foam. 




Fig. 1: Various situations for the boundary of a 2D foam, (a): 
infinite or periodic foam, (b): free cluster (figure taken from 
|25| V (c): cluster of bubbles within a larger foam, (d): foam 
enclosed in a container. Situation (a) is referred to as an un- 
bounded foam, while situations (c) and (d) are referred to 
as bounded foams. Situation (b) is either an unbounded or 
bounded foam, depending on whether or not the surrounding 
air is included as a supplementary bubble. 

Rigorously, the accessible microstates do not all corre- 
spond to the same total surface energy. It is tempting, by 
analogy with the statistical mechanics of a gas, to restrict 
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the equiprobability hypothesis to the microstates of equal 
energy. This refinement appears to be unnecessary: pre- 
vious studies have shown that different arrangements of a 
large number of bubbles of given areas do not affect the 
energy very much 2 : in their computational studies of 2D 
foams under shear, Jiang et al. [27] reported energy fluctu- 
ations less than 2%. Graner et al. [25 obtained the same 
results, both numerically and experimentally: the energy 
values of the different metastable states of a 2D foam lie 
within 2% of the "ground state" value. Kraynik et al. [12] 
conducted similar studies on 3D foams, and reported fluc- 
tuations below 4% for the surface energy of every bubble. 
These observations are also consistent with the so-called 
"equation of state" of a foam [25|l3TH33] , which relates sur- 
face energy E, areas and pressures Pi within a cluster 
of bubbles: E = (Pi — Pq) A^ where Po is the exter- 
nal pressure. In the limit of fixed bubble pressures, the 
total surface energy of such a cluster would be conserved. 
Hence, to a first approximation, one can reasonably as- 
sume that the energy of a large cluster does not depend 
on its configuration, but is directly determined from the 
number of bubbles Nb and the area distribution p(A). 

On a macroscopic scale, the state of the foam should be 
described by a limited number of independent variables 
(besides the number of bubbles and the size distribution). 
These quantities must be "constants of motion" [30], i.e: 
they must keep constant values throughout the dynamics 
of the system (here, a slowly agitated foam). An impor- 
tant observation is that one can define two independent 
quantities which are conserved for an unbounded foam. 
The first one is the total number of sides N s = Xli n i5 
where rii is the number of sides of bubble i. One has 
N s = 27Vp5 for an unbounded foam (7V S , and not Np^, is 
an extensive and fluctuating variable for a finite cluster). 
Thus, according to Eq. (|3]), N s = 6(N B - c). The second 
one is the total cell curvature K tot = J2i K i-> where ki is 
the cell curvature of bubble i, defined as: 

M(i) denoting the neighbouring cells of bubble i. Obvi- 
ously, K to t is additive. Since the terms cancel 
in pairs in the double sum, yielding n to t — for an un- 
bounded foam. The property correlates the 
two adjacent bubbles that share a common edge. It is 
also important to note that - regardless of this mathe- 
matical property - the constraint K to t = is also imposed 
by the curvature sum rule (j2]), as it is illustrated on Fig. 
[2j This rule, which correlates the three adjacent bub- 
bles that share a common vertex, has a physical origin 
(Laplace's law). Although it is unclear whether N s and 
Kt t are the only constants of motion for an unbounded 
foam, the constraints N s = 6(Nb — c) and K to t — must 
be taken into account in the statistical description of such 

2 Actually, there is also some uncertainty on the value of the en- 
ergy of an isolated volume of gas |30| . 



a system. Surprisingly enough, the constraint on the total 
curvature has always been ignored in previous theoretical 
models Q2H22]. 




Fig. 2: Illustration of the equivalence (for an unbounded foam) 
between the sum (over all vertices) S a and the sum (over 
all cells) Ktot = where S a denotes the sum - turning 

clockwise - of the algebraic curvatures Kij + Kjk + km between 
the three bubbles z, j, k that share the same vertex a, while m 
denotes the cell curvature of bubble i (see Eq. Q). Each term 
Kij in the first sum is represented by an arrow pointing from 
i to j. The equivalence between the two sums is immediate: 
the cell curvature m of cell i is represented by the collection 
of arrows coming out of that cell. For an unbounded foam 
at mechanical equilibrium, the curvature sum rule J2]) yields 
S a = at every vertex a, and thus n to t = 0. 

In the light of this discussion, we may restate the pos- 
tulate as: all the accessible states of a given set of Nb 
incompressible cells corresponding to the same values of 
N s and K tot are equally probable, under macroscopic agita- 
tion. By analogy with thermal physics, the microcanonical 
ensemble (Nb, Ktot-, N s ) refers to a large number of copies 
of a foam (with given bubble size distribution p(A)) whose 
parameters Nb-, Ktot-, and N s are fixed. 

Thermodynamic limit. — For a bounded foam, K tot 
reduces to the sum of the side curvatures along its bound- 
ary. Both N s and K tot are fluctuating variables, although 
their values are usually restricted within a certain range. 
For instance, for a free cluster, N s < 2Npb and K to t < 
since the pressure in any bubble of the cluster is higher 
than the external pressure [31]. For a cluster within a 
larger foam, N s < 27Vp5, but n to t can be positive or 
negative. For a foam enclosed in a container, K t ot — 
(the container walls can be regarded as zero curvature 
sides), but N s fluctuates. One can argue whether or 
not the fluctuations of N s and n to t become negligible in 
the thermodynamic limit (Nb — >• oo, total surface area 
A to t —> 00 ? but (A) = A to t/NB remains finite). The av- 
erage number of sides of a large cluster (Nb ^> 1) scales 
as (TVs) °° Nb, while the total curvature, equal to the 
sum of the side curvatures along the cluster boundary, 
scales as (n to t) ~ ^/Nb/(A). The relative standard devi- 
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— 1/2 

ations of such quantities scale as N B 



[30]. Thus, the 



fluctuations of N s and K to t become vanishingly small as 
Nb —> oo. Moreover, their thermodynamic limit values 
are known: (N s ) /N B -» 6, and (fttot) /Nb -> 0. Note 
however that (N*)/N B ~ N B ° converges much faster than 

«ot)/N B 



N B ~ 1/2 . 



Idealized foam. — In order to obtain the probability 
distribution of the number of sides for a cell of given size, 
we need to enumerate the microstates which correspond to 
the same macroscopic state. Ideally, a microstate is spec- 
ified by the curvature, length, position and orientation 
of each Plateau border. However, this description is too 
cumbersome to handle. Moreover, even for a given foam 
topology, the number of different accessible microstates 
depend in a non-trivial way on boundary conditions (e.g.: 
free, periodic, or enclosed cluster) and bubble size distri- 
bution [2 



(a) 




Fig. 3: (a): Real two-dimensional foam under shear, (b): 
Idealization of the foam (mean- field approximation): adjacent 
bubbles i, j and k are now regular cells surrounded by an ef- 
fective foam and disconnected from each other. 



We shall use a simplified microscopic description of the 
foam: we consider that for each particular bubble z, the 
rest of the foam can be replaced on average by a "mean 
field" of identical bubbles. Thus, the bubble i is a regular 
cell with identical curved sides joining at each vertex with 
angles of 120°, as sketched on Fig. [3}}. The cell curvature 
of a n-sided regular cell with area A is [22) [25]: 



6) 



, . v 7r n (n 

n{A,ri) = -— = 

3 e(n)VA 



(5) 



-1/2 



IS 



wWp p(n\ — ^\ n - Q \ ( 1L _ 2L _ sin(7r/n-7r/6) sin(7r/6) \ 

wnereein;- ^ y n 6 S in(7r/n) J 

the elongation of the cell (ratio of perimeter to square- 
root of area). e(n) is a slowly decreasing function, lying 
between e(3) — 3.74 and e(oo) ~ 3.71 [25]. It can be 
noticed that while the (surface) energy of such a bubble 
(~ e(n)\^A) is almost independent of n, its cell curva- 
ture k (A, n) increases rapidly with n, and is not upper- 
bounded. 

In this idealized foam description, each bubble is "dis- 
connected" from the others. By construction, all the phys- 
ical, geometrical and topological constraints are satisfied, 



except the curvature sum rule (|2]) and the Euler-Plateau 
relation (|3]), which both involve adjacency of the bubbles. 
Consistent with the discussion above, the accessible mi- 
crostates of a foam tiling the entire plane are those which 
satisfy N s = 6(7V# — c) and K to t = (the number of bub- 
bles is large enough so N s and K to t can be treated as con- 
tinuous variables). A microstate of the idealized foam is 
specified by the numbers of sides {rii} of all the bubbles. 
Hence, the number of accessible microstates is obtained by 
enumerating the distributions of the N s = 6(Nb — c) sides 
over the Nb bubbles which satisfy the constraint n to t — 0. 
This number is not easy to evaluate, and a grand-canonical 
description shall be more appropriate. 

Grand-canonical ensemble [Nb, , p) . — Con- 
sider then a sample of N B bubbles in the unbounded foam 
(the asterisk denotes the variables of this grand-canonical 
ensemble). This sample can be a cluster of bubbles, or a 
collection of isolated bubbles. Let p*(A) be the distribu- 
tion of bubble size within this sample. Although N B is 
fixed, this system can exchange sides and curvature with 
the rest of the foam, through Tl events. Hence, the to- 
tal number of sides TV* and the total curvature K,* ot are 
now internal variables free to fluctuate for this system 3 . 
We assume that possible other variables required to de- 
scribe the macroscopic state remain fixed. Surface energy 
in particular does not fluctuate, since it depends only on 
N B and the distribution p*(A), and not on the specific 
configuration of the system. Suppose that the rest of 
the foam is large in comparison with the system, so that 
it constitutes a reservoir of sides and curvature. Using 
the formalism of conventional statistical mechanics [30] , it 
comes that the probability for the system to be in the mi- 
crostate (ni, ri2, . . . , rijv* ) is proportional to e~P K tot+v> N s ? 

with K$ ot = Ys?Ji K(Ai, rii) and A^ s * = Y^fji n i- P~ l and V 
(rigorously, /i/3 _1 ) denote respectively the "temperature" 
of the reservoir of curvature, and the "chemical potential" 
of the reservoir of sides. A large number of copies of a foam 
whose parameters (Nb, , /i) are fixed shall be referred 
to as a grand- canonical ensemble. As noticed before, the 
cell curvature of a bubble is not upper-bounded as n — >> oo, 
what ensures that the temperature f3~ l is always positive 
[30] . Finally, the probability for a given cell of size A to 
have n sides is: 



PaH = x _1 (^)e"^ (A ' n)+/in , 



(6) 



where x(A) = E n >3 e~^< A ^ + ^ n denotes the parti- 
tion function of the cell. The average total cell cur- 
vature and average number of sides are, respectively, 
(K* ot ) = -dlnZ/d/3 and (A^ s *) = dlnZ/dp, with InH = 
N B j^°p*(A) lnx(A)dA. E is the partition function of 
the system. Using the formalism of conventional statisti- 
cal mechanics [30] , one concludes that the thermodynamic 



3 iV* and K>t ot are independent variables: there are different ways 
of distributing TV* sides over Ng bubbles; each of these distributions 
corresponds to a different value of K% ot . 
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Table 1: Statistical description of an ideal gas and a 2D foam. 





ideal gas 


Ayj to am 




source of 


Brownian 
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ergodicity 


motion / 


shear /ll 




IT " 

collisions 


events 


constants 
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areas 


degrees of 


positions and 


side numbers 


ireeaom 


momenta 
ir»,p*> 


W 


fixed 


energy E, 


total curvature 


parameters 


number of 


Ktoti number of 


(microcanoni- 


molecules N, 


sides N s , 


cal 


volume V 


number of 


ensemble) 




bubbles N B 


fixed 


temperature T, 


effective 


parameters 


chemical 


temperature 


(grand- 


potential fi, 


effective 


canonical 


volume V 


chemical 


ensemble) 




potential /i, 
number of 
bubbles 



equilibrium of a foam in contact with a reservoir of curva- 
ture and sides coincides with the minimum of the thermo- 
dynamic potential <£> = — /3 _1 ln£. The analogy between 
the statistical descriptions of an ideal gas and a 2D foam 
is summarized in Table HJ 

As noted by Graner and coworkers [25], the elongation 
of a regular cell is almost constant: e(n) 3.72. With 

this simplification, the total surface energy of the cluster 
is strictly conserved (i.e., it does not depend on the cluster 
configuration). Besides, the distribution © simplifies to 

p A (n) = X '- 1 (A)e-^-^'^- 6 \ (7) 

with X ' = Xe" 6/i , P'{A) = 7r(3/(3eVA) and fi'(A) = 
\± — 2ir f3 1 '(ey/A) . This is the exact distribution intuited 
independently by Schliecker and Klapp [T7J[T8] and Sher- 
rington and coworkers 281129]. except that here j3' and p! 
explicitly depend on the bubble size A. Such a dependence 
is necessary to reflect the correlations between bubble size 
and bubble shape which have been observed experimen- 
tally [5,6 . It can be noted that the average number of 
sides n(A) = X^n>3 n ^( n ) mcreases with A. This re- 
sult is consistent with experimental observations [Tjl39]: 
larger bubbles have more sides since they are surrounded 
by smaller bubbles. 

It must be pointed out that the present theory contains 
no free parameters in the thermodynamic limit: (N*) /Ng 
and (f^t ot ) /N B have known values, and f3 and \i are di- 
rectly obtained by extremizing the grand-canonical en- 
tropy per bubble 

S(/3,/x) = r P*(A)\n X (A)dA + f3 { -^ - (8) 

with (N;) /N* b = 6 and ( K $ ot ) /N* B = 0. 



Discussion. — The expression ([7]) allows to de- 
duce the distribution of number of sides per cell p(n) 
from the distribution of bubble size p(A) [p(n) = 
p(A)pA(n)dA], and thus to correlate topological and 
geometrical disorders. Let us study the implications 
of the present theory with the simple case of an infi- 
nite/periodic monodisperse foam: p(A) = 5 (A — Ao). In 
that case, the thermodynamic limit values (N*)/Ng = 6 
and (K% ot )/Ng = give, respectively, (n — 6) =0 and 
((n — 6) 2 ) = 0. Thus, the distribution of side number tends 
to the Kronecker delta distribution: p(n) = PA (n) = S n ^. 
As expected, all the bubbles of an unbounded monodis- 
perse foam have hexagonal shape. 
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Fig. 4: Theoretical and experimental distributions p(n). The 
experimental distribution is taken from 5 . The theoreti- 
cal distribution is obtained by fitting the function p(n) = 
J °° p(A)pA(n)dA to the data. Inset: topological disorder 
An/(n) vs geometrical disorder A A/ (A) for four different 
foams. Black dots: obtained from theory; red curve: power- 
law fit a (AA/{A)) b , with a = 0.30 ± 0.03 and b = 0.74 ± 0.08; 
blue curve: linear fit a A A/ (A) + 6, with a = 0.30 ± 0.02 and 
b = 0.03 ±0.01. 

We also compare the model with the four experimen- 
tal distributions reported in the Quilliet et at. paper [5]. 
Foam samples contain a few hundreds of bubbles. This 
number is large enough to consider that (N*)/Ng ~ 6, 
but not large enough to assume that (K*ot)v (A*) /N B — 0- 
Since the finite value of (f^tot)/^B * s unknown, j3 and \i 
cannot be obtained by extremizing S{fi,\L). Instead, they 
are obtained by fitting the theoretical distribution of sides 
p(n) = J °° p(A)pA(n)dA - in which p(A) is the experi- 
mental size distribution, and pA(n) is given by Eq. ([7j) 
- to the data. Regression is performed under the con- 
straint Xln>3 n P( n ) — 6, so there is only one adjustable 
parameter really. Fig. H] compares a theoretical distribu- 
tion p(n) obtained this way with the corresponding exper- 
imental distribution taken from [5] . The theoretical curve 
reproduces the experimental data well. Comparison of the 
model with the other distributions of [5] (not shown here) 
gives similar results. We also plotted the relative standard 
deviation An/(n) = yj (n 2 ) — (n) 2 /(n) of the theoretical 
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distributions, as a function of the relative standard devia- 
tion A A I (A) of the size distributions for the four samples 
(inset of Fig. 2]). A power-law fit of this plot is in good 
agreement with the relationship found by Quillet et al: 
An/(n) = 0.27(AA/(A)y-*. However, it can be noticed 
that the power-law fit can be hardly distinguished from a 
linear fit. 

In summary, we developed a theoretical model to de- 
scribe the state of a two-dimensional foam under slow ag- 
itation, using a formulation closer to conventional statis- 
tical mechanics than information theory. We show that 
the total number of sides and the total cell curvature 
- rather than energy - are the relevant variables to de- 
scribe the macroscopic state of a foam. The distribution 
of sides of a cell of given size is derived. This result allows 
to correlate the size and shape distributions. Theoreti- 
cal size-topology relations deduced from the theory are 
in very good agreement with existing experimental data. 
However, a free parameter has been added for the com- 
parison with experiments, due to the indeterminacy of 
the experimental values of (f^tot)/^B- Furthermore, the 
grand-canonical description requires that 1 <C Ng <C Nb- 
This condition cannot be checked on the available data. 
Further experiments taking these considerations into ac- 
count should allow to confirm the validity of the present 
(zero-free-parameter) model. 

The formalism developed here can be extended to three- 
dimensional foams, and to coarsening (and coalescing) 
foams. Coarsening has consequences for the size- topology 
correlations particularly, since the average area of n-sided 
cells increases more rapidly in that case [6,39,40]. Exten- 
sion of this formalism to other cellular systems (concen- 
trated emulsions, cell aggregates) is also under investiga- 
tion. 



I wish to thank J. B. Fournier, F. Graner and F. Van 
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